Higher order finite difference schemes for the magnetic induction equations

We describe high order accurate and stable finite difference schemes for the initial-boundary value problem associated with the magnetic induction equations. These equations model the evolution of a magnetic field due to a given velocity field. The finite difference schemes are based on Summation by Parts (SBP) operators for spatial derivatives and a Simultaneous Approximation Term (SAT) technique for imposing boundary conditions. We present various numerical experiments that demonstrate both the stability as well as high order of accuracy of the schemes.Comment: 20 page

Galaxies in box: A simulated view of the interstellar medium

We review progress in the development of physically realistic three dimensional simulated models of the galaxy.We consider the scales from star forming molecular clouds to the full spiral disc. Models are computed using hydrodynamic (HD) or magnetohydrodynamic (MHD) equations and may include cosmic ray or tracer particles. The range of dynamical scales between the full galaxy structure and the turbulent scales of supernova (SN) explosions and even cloud collapse to form stars, make it impossible with current computing tools and resources to resolve all of these in one model. We therefore consider a hierarchy of models and how they can be related to enhance our understanding of the complete galaxy.Comment: Chapter in Large Scale Magnetic Fields in the Univers

Diffusion in supersonic, turbulent, compressible flows

We investigate diffusion in supersonic, turbulent, compressible flows. Supersonic turbulence can be characterized as network of interacting shocks. We consider flows with different rms Mach numbers and where energy necessary to maintain dynamical equilibrium is inserted at different spatial scales. We find that turbulent transport exhibits super-diffusive behavior due to induced bulk motions. In a comoving reference frame, however, diffusion behaves normal and can be described by mixing length theory extended into the supersonic regime.Comment: 11 pages, incl. 5 figures, accepted for publication in Physical Review E (a high-resolution version is available at http://www.aip.de./~ralf/Publications/p21.abstract.html

Simulation techniques for cosmological simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.Comment: 42 pages, 16 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 12; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke

Magnetic Field Amplification in Galaxy Clusters and its Simulation

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.Comment: 60 pages, 19 Figure

Coronal Magnetic Field Evolution from 1996 to 2012: Continuous Non-potential Simulations

Coupled flux transport and magneto-frictional simulations are extended to simulate the continuous magnetic-field evolution in the global solar corona for over 15 years, from the start of Solar Cycle 23 in 1996. By simplifying the dynamics, our model follows the build-up and transport of electric currents and free magnetic energy in the corona, offering an insight into the magnetic structure and topology that extrapolation-based models cannot. To enable these extended simulations, we have implemented a more efficient numerical grid, and have carefully calibrated the surface flux-transport model to reproduce the observed large-scale photospheric radial magnetic field, using emerging active regions determined from observed line-of-sight magnetograms. This calibration is described in some detail. In agreement with previous authors, we find that the standard flux-transport model is insufficient to simultaneously reproduce the observed polar fields and butterfly diagram during Cycle 23, and that additional effects must be added. For the best-fit model, we use automated techniques to detect the latitude–time profile of flux ropes and their ejections over the full solar cycle. Overall, flux ropes are more prevalent outside of active latitudes but those at active latitudes are more frequently ejected. Future possibilities for space-weather prediction with this approach are briefly assessed

Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics

This article serves as a summary outlining the mathematical entropy analysis of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD equations as they are particularly useful for mathematically modeling a wide variety of magnetized fluids. In order to be self-contained we first motivate the physical properties of a magnetic fluid and how it should behave under the laws of thermodynamics. Next, we introduce a mathematical model built from hyperbolic partial differential equations (PDEs) that translate physical laws into mathematical equations. After an overview of the continuous analysis, we thoroughly describe the derivation of a numerical approximation of the ideal MHD system that remains consistent to the continuous thermodynamic principles. The derivation of the method and the theorems contained within serve as the bulk of the review article. We demonstrate that the derived numerical approximation retains the correct entropic properties of the continuous model and show its applicability to a variety of standard numerical test cases for MHD schemes. We close with our conclusions and a brief discussion on future work in the area of entropy consistent numerical methods and the modeling of plasmas